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ABSTRACT 


The  problem  of  an  axial  slot  on  a  circular  cylinder  clad  with  a 
radially  inhomogeneous  plasma  (dielectric)  was  considered  recently  by 
Rusch  [l964] .  Although  his  formulation  is  quite  general  inasmuch  as  it 
admits  a  realistic  approximation  to  the  plasma  sheath  inhoraogeneities,  the 
resulting  solution  is  not  well  suited  for  computations  of  field  patterns 
of  large  cylinders. 

In  the  present  formulation  we  consider  the  dielectric  permittivity 
profile  of  the  form  e(p)  =  (p/b)^,  a  <  p  <  b  ,  where  'a'  is  the  radius  of 
the  conducting  cylinder,  *b'  is  the  outer  radius  of  the  dielectric  coating 
and  'p'  is  an  arbitrary  parameter.  It  follows  that  the  assumed  profile  is 
capable  of  representing  an  inhomogeneous  cold  plasma  when  p  >  0  and  a  real 
dielectric  when  p  <  0.  The  apparent  advantage  of  the  present  formulating 
is  the  fact  that  the  wave  equation  can  be  solved  in  terms  of  known  functions 
and  the  solution  can  be  extended  to  large  cylinders. 

Field  expressions  appropriate  to  small  and  large  cylinders  are 
found  using  standard  methods  of  harmonic  series  representation,  Watson  Trans¬ 
formation  and  saddle  point  integration.  In  the  case  of  large  cylinders  coated 
with  cold  plasma,  the  radiation  patterns  are  plotted  for  various  combinations 
of  the  cylinder  radius,  thickness  of  the  coating,  and  the  inhoraogeneity  gra¬ 
dient.  The  radiation  patterns  are  in  good  agreement  with  qualitative  argu¬ 
ments  of  geometric  optics. 
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1.  Introduction 


The  problem  of  a  radiating  slot  on  a  dielectric  clad  cylinder  occu¬ 
pies  a  relatively  important  position  in  microwave  engineering.  Its  importance 
is  due  primarily  to  the  phenomenon  of  a  plasma  sheath  formed  around  a  hyper¬ 
sonic  vehicle  during  its  reentry  phase.  Since  the  plasma  can  be  represented 
under  certain  conditions  by  it a  equivalent  dielectric  permittivity,  the  problem 
of  a  plasma  covered  antenna  can  be  reduced  to  one  of  a  dielectric  clad  antenna 
that  can  be  treated  as  a  well-defined  boundary  value  problem. 

The  boundary  value  problem  approach  to  the  plasma  sheath  phenomenon 
was  used  in  the  past  by  several  authors.  Tamir  and  Oliner  [l962]  and  Omura 
[l962]  obtained  theoretical  radiation  patterns  of  a  dielectric  clad  slot  on  ? 
ground  plane.  Tyras  et  al  [l965]  obtained  experimental  confirmation  of  their 
results  by  means  of  an  artificial  plasma  simulation  technique.  The  problem 
of  an  axial  slot  on  a  circular  cylinder  clad  with  a  dielectric  layer  was  con¬ 
sidered  by  Hasserjian  [l965] .  He  formulated  the  solution  in  terms  of  an 
integral  representation  from  which  he  was  able  to  find  the  far  field  pertinent 
to  a  cylinder  of  a  large  radius. 

The  problem  of  a  slot  clad  with  inhomogeneous  dielectric  layer  was 
attacked  by  two  different  methods.  In  the  first  method,  the  inhomogeneity  of 
the  layer  was  approximated  by  a  series  of  discrete  homogeneous  and  parallel 
layers  [Hamm  and  Tyras,  1965;  Harris  and  Pachares,  1965;  Harris,  1965].  The 
resulting  formulas  are  necessarily  long  and  complicated  but  they  can  be  pro¬ 
grammed  for  a  digital  computer  to  obtain  numerical  results.  In  the  case  of 
cylindrical  geometry,  the  formulas  obtained  by  this  method  converge  slowly 
when  the  radius  of  the  cylinder  is  large  and,  consequently,  their  application 
is  limited. 
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Another  method  of  approach  to  the  problem  of  inhomogene i ties  in  the 
layer  is  that  of  a  rigorous  solution.  Inasmuch  as  no  general  solution  to  the 
differential  equation  describing  wave  propagation  in  inhomogeneous  regions 
exists,  this  approach  suffers  from  the  lack  of  versatility.  However,  Rusch 
[l964]  found  a  solution  to  an  inhomogeneous  profile  that  allows,  perhaps,  the 
most  realistic  approximation  to  an  inhomogeneous,  continuously  radially 
stratified  plasma  sheath.  He  chose  a  model  of  the  dielectric  permittivity 
profile  of  the  form  e(p)  =  Ap  +  Bp  +  C  where  A,  B,  and  C  are  arbitrary 
constant.  The  two  linearly  independent  solutions  of  the  wave  equation  are 
constructed  using  the  well-known  method  of  Frobenius.  Consequently,  his 
final  expression  for  the  fields  appears  as  an  infinite  series  that  cannot  be 
expressed  in  terms  of  known  functions.  When  the  radius  of  the  cylinder  is 
small,  the  resulting  series  converges  well  and  numerical  results  can  be  ob¬ 
tained.  When  the  radius  of  the  cylinder  is  large,  which  is  the  case  of  great 
practical  importance,  the  resulting  series  is  slowly  convergent  and  it  is  not 
suitable  for  numerical  calculations. 

In  this  paper  we  choose  the  dielectric  permittivity  profile  of  the 
form  e(p)  =  (p/b)^  where  * b *  is  the  outer  radius  of  the  coating  and  ‘p*  is 
an  arbitrary  constant.  This  partJcular  model  admits  exact  solution  of  the 
pertinent  differential  equation  in  terms  of  Bessel  functions.  The  model  is 
capable  of  representing  an  increasing  or  decreasing  profile  when  * p *  is 
positive  or  negative  respectively.  When  *p*  is  positive,  this  model  can 
be  related  to  the  plasma  frequency  in  the  form  (w0/w)^  «  1  -  (p/b)^  . 

It  will  follow  that  this  representation  will  admit  the  specification  of  the 
plasma  frequency,  to0,  between  the  limits  of  the  operating  frequency  and 


zero. 
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The  principal  advantage  of  the  present  formulation  is  the  fact  that 
the  solution  can  be  expressed  in  terms  of  known  functions.  After  obtaining 
the  familiar  harmonic  series  representation,  the  integral  representation  can 
be  found  using  the  well-known  Watson  Transformation.  The  latter,  in  turn, 
enables  obtaining  the  far  field  expression  for  an  arbitrarily  large  cylinder 
which  seems  to  be  of  some  practical  importance. 

2.  Theoretical  Development 

Consider  an  axially  slotted  cylinder  clad  with  an  concentric  in¬ 
homogeneous  layer  of  a  dielectric  that  is  a  function  of  the  radial  variable 
only.  The  geometry  of  the  problem  is  shown  in  Fig.  1.  In  the  inhomogen-.  us 
dielectric  region  the  fields  must  satisfy 

**>■  -f 

Vx  E|  °  io^0Hl 

Vx  Hi  «>  iw€0€(p)  Ex 

E$l(a,*)  a  1  &(*)  (1) 

**  i  t 

where  the  harmonic  time  dependence  of  the  form  e  has  been  assumed  and 
suppressed  throughout.  The  excitation  considered  is  of  the  form  of  an  in¬ 
finitely  thin  slot  fed  by  a  voltage  V  (delta  function  slot).  In  the  free 
space  region  the  fields  must  satisfy 

V  x  E0  =  iw  p0  hJ) 

V  x  Hq  **  -iw€0^,  .  (2) 

Moreover,  at  the  interface  p  =  b  the  tangential  components  of  the  electric 
and  the  magnetic  fields  must  be  continuous  which  implies 
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H*1  "  Kzo 

^❖1  =  >  (3) 


since  the  assumed  geometry  does  not  admit  any  other  tangential  field  components. 

Performing  a  curl  operation  on  the  second  equation  of  (1)  and  sub¬ 
stituting  the  first,  one  obtains  a  differential  equation  in  terms  of  the 
magnetic  field  only,  that  is 


V2  5^  +  i  7  €x  4hl  +  €  Hj.  o  0 


(4) 


where  we  used  a  well-known  vector  relation.  The  axial  component  of  the  vector 
equation  (4)  can  be  shown  to  be 


i  a  aHzU  i  a2Hzl  i,  de  at^  i  ae  aBzl  ]2  .. 

o  Sp^p  -§r>  +?'a7“'7(  5Sar_  +  7a7  5T-  >  +  k°e  ^  ■ 0 


(5) 

which  in  the  case  that  the  dielectric  permittivity  is  a  function  of  the  radial 
variable  only  becomes 


1  5  ^^zl\  1  d€  ^Hzi  i  2 

n  aS(o  §T->  "  7  Tp  ST"  +  ~T.  ~ '5Zr  +  k°£Hzi  =  0 


(6) 


Similarly,  it  can  be  shown  that  in  the  free  space  region  the  magnetic  field 
must  satisfy 


ST0 


C>2hZo  . 
c)$2 


k2H 

KOnZO 


(7) 


It  will  be  appropriate  to  seek  the  solution  of  the  form 
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H2(p,$)  =  ^m(p)eim0  (8) 

-co 

from  which  it  will  follow  that  the  radial  functions  ^rai(c)  snd  i^Cc)  mus«- 
satisfy  the  ordinary  differential  equations 


where  tie  denoted  differentiations  are  with  respect  to  p  . 

The  further  success  in  proceeding  with  the  solution  is  contingent 
upon  being  able  to  solve  the  differential  equation  (9).  Since  a  general 
solution  for  an  arbitrary  function  e(p)  is  not  known  to  exist,  one  must  at 
this  point  specify  the  functional  dependence  of  the  dielectric  permittivity, 
e  .  Rusch  [l964 ]  chose  a  model  of  the  dielectric  permittivity  of  the  form 
€ ( p)  =  Ap  r  Bp  +  C  where  A,  B  and  C  are  constants  and  found  two  linearly 
independent  power  series  solutions  to  (9)  using  the  method  of  Frobenius 
[Dettman,  1962],  His  model  is  quite  versatile  inasmuch  as  it  allows  fairly 
realistic  approximations  to  the  electron  density  profile  in  a  plasma  sheath. 
The  price  paid  for  this  is  the  complexity  of  the  final  solution  which  cannot 
be  expressed  in  terms  of  known  functions. 

In  this  paper  we  shall  choose  the  profile  of  the  dielectric  per¬ 
mittivity  in  the  form 

€(P)  =  (§)2P  (11) 


i 
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where  ’bf  is  the  outer  radius  of  the  coating  and  ’p*  is  an  arbitrary  number. 
In  what  follows  we  shall  show  that  this  particular  choice  of  e(p)  results  in 
solutions  to  (9)  in  terms  of  known  and  tabulated  functions. 

The  particular  model  of  e(p)  given  in  (11)  is  capable  of  represent¬ 
ing  and  increasing  profile  when  1 p *  is  positive  and  a  decreasing  profile  when 
*p*  is  negative.  When  'p*  is  positive,  this  permittivity  profile  can  be 
easily  related  to  the  electron  density  profile  in  a  lossless,  cold  plasma  via 
the  relation 


2  7  2 

where  w0  is  th^  plasma  frequency  and  (w0/u>)  =  8.06  x  10  N/f  ,  where  N  is 

the  number  of  free  electrons  per  cubic  centimeter.  The  possible  profiles  of 
(u>0/w)  versus  p/b  are  shown  in  Fig.  2a.  It  will  follow  that  the  present 
representation  will  admit  the  specification  of  the  plasma  frequency  u>0  be¬ 
tween  the  limits  of  the  operating  frequency  w  and  zero.  Once  the  plasma 
frequency  is  specified  at  the  skin  of  the  cylinder,  however,  the  details  of 
the  profile  are  fixed.  The  latter  restriction  is  not  a  serious  one,  however, 
inasmuch  as  it  is  known  that  the  main  contributions  to  the  propagation 
characteristics  will  come  from  the  bulk  effect  of  the  inhomogeneity  rather 
than  from  the  details  of  the  permittivity  profile  [Liu  and  Wetzel,  1965]. 

When  *p'  is  negative  then  the  assumed  functional  dependence  for 
e(p)  in  (11)  will  represent  a  profile  that  decreases  with  an  increasing 
radial  distance  as  shown  in  Fig.  2b.  This  case  is  then  applicable  to  the 
problem  of  a  cylinder  clad  with  real  inhomogeneous  dielectric. 

When  the  form  of  e(p)  given  in  (11)  is  substituted  into  (9),  one 
obtains  the  differential  equation 


-1 
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«Cl  +  +  (ko  ’>"2P  02P  -  3>*m  1  =  0 


(12) 


which  can  be  related  to  the  known  equation  [dver,  1964] 


W"  +  (i^)w‘  +  (\2q2z 


2q-2  +  £?=^ai)H  „  o 


(13) 


having  solutions  of  the  form 

W  =  zP  £v(\zq)  (14) 

where  CyO  is  any  Bessel  function.  Comparing  the  coefficients  of  (12) 
and  (13)  we  find  that  the  solutions  to  (12)  must  have  the  form 

tfrnl  :  PP  ^(Xp^1)  (15a) 

where 


(p+DbP 


v 


sU 


+  p‘ 


p+1 


(13b) 


(15c) 


The  solution  (15a)  is  evidently  valid  for  all  p  ^  -1  and  when  p  =  0  it 
becomes  identical  to  the  corresponding  solution  in  free  space. 


2.1  The  Case  p  ^  -1 

The  complete  solution  can  in  this  case  be  written  in  the  following 


form 
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00 

“  PP  Y  ^  [cI(Jv(XpP+1)  +  DmH^1)(\pPf1)]  (16a) 

•  00 


CO 

Hzo  -  Y  *lm4  F^^OcoP)  .  (16b) 

-  00 


Noting  that  the  tangential  electric  field  components  are  expressed 


as  follows 


E^i  = 


-i  /b\2P  SHzl 
gj€0  (  a)  c)p 


-j  ^^20 
<>0  “  0)€o  dp 


(17a) 

(17b) 


one  obtains,  applying  the  boundary  conditions  at  p  *=  a  and  p  =  b,  the  set  of 
simultaneous  algebraic  equations 


/all 

a12 

l°m 

iwe0V 

I M 
0 

a21 

a22 

a23 

*  2nab2P 

a31 

a32 

a33  j 

°! 

where 

k  p+1 

all  "  b?  1JV(A)  +  ^(“)  Jv(A)] 

a21  “  bP  VB> 

a22  »  bP  H^(B) 

a23  "  -  ^(Bo) 

a31  -  bP[j(,(B)  +  §~  JV(B)] 

a32  a  bP[Hj1)f(B)  +|-  ^(B)] 

a33 


(19) 
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and 


~  ^oa 
Bq  =  kpb 


A 


p  +  1 


B  = 


Bp 

p  +  1 


(20) 


The  system  of  simultaneous  equations  can  be  solved  readily  using 
Cramer's  rule.  The  field  in  the  air  is  found  to  be 


00 

~^£0(p~bl)V  \  -(mS  (kpp) 

■  *2^  6  ^ 


(21) 


where 


-  <JvHll)  '  <Bo) 


(22) 


and 


Pv  -  fJi<A>  +  ^  JV<A>]  fH$1> '  (B)  +  f-  hJ^CB)] 

-  [h^1)  *  (A)  +£-  f-]  41}(A)][j'(B)  +£-  Jv(B)]  (23a) 

B0  \al  Bo 

Qv  -  hv°  <«  fj;<A> 

-  JV(B)  [H$1) '  (A)  +  2-  /-)  ?+1  Hy1^  (A) ]  . 

Bo  \aI 

As  a  partial  check  on  the  above  result,  we  set  p  *»  0  which 


corresponds  to  a  cylinder  in  the  air.  In  this  case  one  finds 
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A  2i  (1)’ 

=  rib  ^  <Ao> 


(24) 


which  upon  substitution  into  (21)  gives 


„  IV  \  iraO  *4^ (koP> 
zo  "  2nZ0a  :  e  5J1) ' (Aq) 


(25) 


-CO 


,  1/2 


where  ZQ  =»  (p0/e Q)  '  is  the  free  space  impedance.  The  expression  in  (25) 
can  be  recognized  as  the  one  appropriate  to  an  axially  slotted  cylinder  in 
free  space  [wait,  1959] . 

Another  partial  check  can  be  performed  setting  b  =  a.  In  this 
case  one  finds 


4n  "  ^ r 1  '  (Ao> 


(26) 


which  when  substituted  into  (21)  again  gives  (25)  as  it  should. 


2.2  The  Case  p  =  -1 

In  the  previous  case  we  have  excluded  the  case  p  =»  -1  because  the 
solution  to  the  radial  wave  function  was  not  defined  therein.  We  shall  now 
examine  this  particular  case  in  detail. 

Putting  p  »  -1  in  (12)  one  obtains 


*2  2 

->  ,  Dn  ~  m 

\^"1  +  1.  +  - - —  'll 


’'ml  1  p  ^ml 


-  V  s  a  0 

2  ml 


(27) 


This  differential  equation  has  simple  solutions 


-1+q  -1-q 

♦ml  :  0  '  p 


1  +  m2  -  Bo2  . 


(28) 
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The  field  in  the  air  now  becomes 

iw€0bV  y-1  eim<>  q(ab>q  (kQp) 

^  112  _4J  '  (Bo)  [(q-Da2q  +  (q+l)b2q]+(l-q2)H^1)  (B0)  [b2c*-a2q] 

(29) 


As  a  partial  check  on  this  result  one  can  set  a  =  b  in  which  case  one  again 
obtains  (25). 

Inasmuch  as  this  solution  is  of  rather  restrictive  application,  it 
will  not  be  pursued  any  further. 

3.  Application  to  Small  Cylinders 
The  formulas  (21)  and  (29)  constitute  the  formal  solution  to  the 
problem  for  any  inhomogeneous  layer  described  by  (11)  and  any  radius  of  the 
cylinder.  In  the  case  when  the  radii  of  the  cylinder  and  the  coating  are 
small,  considerable  simplifications  to  these  formulas  are  possible. 

It  will  be  convenient  to  rewrite  (21)  in  the  form 

C° 

-ue0(p+l)V  —  H<l)(k0p) 

^  -  tha  L  £" cos  — — 

o 

fl  ,  ra  =  o 

era  =  1 

[2  ,  m  t  0  .  (30) 

Now  if  the  radius  of  the  coating  is  sufficiently  small  so  that 

B0  «  P  +  1  (31) 


the  following  approximations  to  the  various  Bessel  functions  are  possible 


12 


1 

r(v+i) 


Jy(X)  ~  £  Jy^x) 


(x) 


-.  -i  r  (v) 

It 


(x)  -  ZZ  hJ1}  (x) 

X 


(32) 


Using  these  approximations  in  (21),  one  obtains  its  approximate 
forn  valid  for  B0  «  p  +  1,  p  ^  -1  . 


H 


u>e0V 


ZO  4  *>  o 


-  B(f)pV  (If  C" 
1 


m 


(33) 


where 


Rju  =  (tn-p-\'ra2+p2)  (|)^m  +P#"  -  (m-p  -r  s/m^+p^)  (|)  .  (34) 


It  will  be  noted  that  the  first  term  in  (33)  correctly  represents 
the  field  of  an  isolated  magnetic  current  line  source  of  the  strength  I^«  -V 
The  remainder  of  the  terras  in  (33)  evidently  represents  the  correction  terms 
for  the  non-zero  radius  of  the  cylinder  and  the  effect  of  the  coating. 

4.  Application  to  Large  Cylinders 
When  the  radii  of  the  cylinder  and  the  coating  are  not  small,  the 
harmonic  series  representative  of  the  solution  derived  earlier  is  slowly 
convergent  and  consequently  many  terms  are  needed  for  a  given  degree  of 
accuracy.  In  this  case,  an  alternate,  integral  representation  of  the  solu¬ 
tion  can  be  found  using  the  well-known  Watson  Transformation.  Subsequently, 
this  integral  representation  can  be  used  to  find  the  "residue  series  repre¬ 
sentation"  valid  in  the  shadow  region  and  the  "geometric  optics  field"  valid 


in  the  illuminated  region. 


4.1  Integral  Representation  of  the  Solution 


If  the  denominator  AJ(  in  (21)  has  no  zeroes  on  the  real  axis  in  the 
complex  p-plane,  then  it  will  follow  that  the  harmonic  series  representation 
(21)  can  be  put  in  the  form 


H 


zo 


mc0(p-fl)V  PH^Ocpp) 

12 AqBq  J  A. 

c 


eip(C-n) 

— - -  dp 

sinpn 


(35) 


where  C  is  a  contour  in  the  complex  p  -  plane  that  runs  slightly  above  and 
slightly  below  the  re-1,  axis  in  the  clockwise  direction  in  such  a  way  that 
it  encloses  all  of  the  zeroes  of  sin  M  n  .  Alternately  (35)  can  be  written 
in  the  form 


^zo 


m 

m€q(p+1)V  r  H^i;(kQp) 

ijt^AoB0  J  A^_ 

-00 


cos  p($-Jt) 

— - -  dp 

sm  p  n 


(36) 


where  the  path  of  integration  is  understood  to  run  slightly  above  the  real 
axis . 


4.2  Fields  in  the  Shadow  Region 

In  order  to  find  the  residue  series  representation  we  must  evaluate 
the  integral  in  (36)  by  closing  the  contour  with  a  large  semi-circle  in  the 
upper  half  plane  enclosing  all  of  the  poles  produced  by  the  zeroes  of  A^. 

To  this  end  we  must  investigate  the  behavior  of  the  integrand  in  (36)  for 
large  values  of  p. 

Let  us  devide  the  upper  half-plane  as  follows: 

Region  I:  0  <  arg  p  <  n/2  -  5 

Region  II:  n/2  -  5  <  arg  p  <  n/2  +  & 

Region  III:  n/2  +  6  <  arg  p  <  n 
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and  note  the  following  asymptotic  behavior  of  the  Bessel  functions  [Honl  et  al, 


196 1] 

I: 

h£3  )  (x)~ 

Y  «P  l ex/ 

(«) 

(37a) 

II: 

h£L)  (x)~ 

A m'v  , 

I  jtp  lex/ 

(37b) 

III: 


h£L)(x)~ 


Ju(x)~ 

V-2np 


(37c) 


Using  these  asymptotic  forms,  (22)  can  be  written 


<Vi 


(Am)h 


12u  /bW 
nAB0  ( a / 


i2i 


nAB, 


-p 


ivP2  x 
2AB0m  av 


H^1} (Bc) 

hJ1} (B0) 

H(1)(Bn) 

p  o 


where 


Xv  -  Jv(A)Yv(B)  -  JV(B)YV(A). 


(38a) 


(38b) 


(38c) 


(39) 


The  asymptotic  evaluation  of  Xv  in  (39)  is  somewhat  cumbersome 
inasmuch  as  the  straightforward  application  of  (37c)  yields  a  value  of  zero. 
In  order  to  obtain  a  more  accurate  estimate,  one  may  apply  the  uAiltiplica- 
fion  theorem  for  the  Bessel  functions  [olver,  1962] 

00 

tv(A)  ■  ev(OB>  =«V^I~r-(f)n  tv+nfW 

n=o 


(40a) 


with 


“  “  ft) 


(40b) 


Substituting  (40a)  into  (39)  gives 


Xy  =  a 


in  / B\n 


n!  (2 


)  fYv(B)  Jv+JB)  -Jv(B)Yv4fl(B)] 


Now  using  the  Wronskian 


VB>WE)  -  Vl<»V‘)  *  ni 


and  the  recurrence  relation 


ew2<B>  -  fi  (B>  ■  <VB) 


which  for  a  large  v  can  be  approximated  by 


~(r)  U i<« -(f) ’cvcb), 


one  can  sura  the  series  in  (41)  to  obtain 

V  _  f„v(1-a2)  n 

xv  ~  m  <e  '  ' 


In  region  III,  p  and  consequently  V  will  have  a  negative  real  part 
so  that  the  exponential  in  (45)  will  have  negligible  contribution.  Thus 


(38c)  becomes 


(a  ) - =*■ 

'/VlII  2nAB, 


IP2  la\V  M) 


AD 
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Setting  p 

that  for  a  large  R 


Hi1)(koP)c°s  »*(*-it)dn 


sin  p  n 


R  exp  (i©)  on  the  semicircle,  it  is  now  easy  to  show 

^  e-R[ysin0  +  cos©  in  (p/a)] 

-Rfysin©  -  cos©  in  (p/a)1 

^  6 

^  R2e-R[ysin0  -  cos©  in  (p/a)]  (47) 


in  regions  I,  II,  and  III  respectively,  where  y  =  o  for  0  <  <>  <  «  and 
y  *=  2n  -  0  for  n  <  $  <  2it  .  Thus  the  contribution  from  the  large  semicircular 
arc  vanishes  everywhere  in  the  upper  half  plane. 

The  integral  in  (36)  can  now  be  evaluated  to  give 


_  2ue0(ptI)V  Y 

Hzo  =  *A0B0  /  Cn 

n=l 


u(1)n,  ncos 
Hpn  (k0P>- 


sxn  pnn 


(48a) 


where 


(48b) 


and  lu  are  the  zeroes  of  A  =0. 
n  y  n 

In  the  far  zone,  the  Hankel  function  can  be  replaced  by  its 
asymptotic  form 


H 


(1) 

Pn 


ei(kcp  -  jt/4 


pnrt/2) 


(49) 


and  noting  that 


cos  Mn(»“rt)  e-iMnn/2  ^ 
sin  pnn 


-i  [eiMn<3jt/2"*>+e  W*'"/2)] 


(50) 
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it  will  follow  that  the  series,  in  (48)  will  converge  rapidly  in  the  shadow 
region  of  the  slot,  it/2  <  «  <  3it/2. 

The  zeroes  of  Apn  =  0  are  evaluated  approximately  in  the  Appendix. 

4.3  Fields  in  the  Illuminated  Region 
To  find  the  fields  in  the  region  directly  illuminated  by  the  slot 
we  must  evaluate  the  integral  in  (36)  at  its  stationary  points.  In  what 
follows  it  will  be  shown  that  these  stationary  points  will  be  found  on  the 
real  axis.  In  order  to  assure  freedom  in  deforming  the  contour  of  integra¬ 
tion  so  as  not  to  interfere  with  the  poles  of  the  integrand  due  to  the  zeroes 
of  sin  p  it,  it  will  be  convenient  to  employ  the  transformation 

cos  p(»-it)  e  co_s  Hi  gipit  .  ieiMe  (51) 

sin  p  n  sin  pit 


which  enables  writing  (36)  in  the  form 


H 


zo 


-oie0(p+l)V 

jt^AoBq 


(ill  +  l2> 


(52) 


where 


_ck/ 


(kpP> 

M 


cos  E  1  ipit 
sin  p  it  e 


dp 


(53) 


and 


^p 


e1^ 


dp  . 


(54) 


The  integral  1^  can  be  readily  evaluated  by  closing  the  contour 
with  a  large  semicircle  in  the  upper  half-plane,  the  contribution  of  which 
is  negligible,  and  integrating  about  the  poles  produced  by  the  zeroes  of 
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A,,  =0.  The  result  is 


*1 


nVt 


cos  i~n? 
sin  pnit 


e 


i^n^ 


(55) 


Using  the  asymptotic  form  of  the  Hankel  function  in  (49)  and  a  simple  trans¬ 
formation,  one  can  readily  show  that  in  the  far  2one 


h 


i(k0p-n/4)y 

n=l 


r  iun(3«/24-5)  i^n(3it/2-4>) 

e. _ -re  _ 

1  -  ei2jUin 


(56) 


Since  pn  has  a  positive  imaginary  part  it  will  follow  that  the  contribution 
of  1-^  in  the  illuminated  region  will  be  small. 

To  evaluate  the  integral  l£  we  first  note  the  following  asym;  cotic 
forms  valid  for  a  large  value  of  the  argument  [Konl  et  al,  19 6 1] 


h^cv- 


/I  (b^v  1/4  |M|  <  bo 


(57a) 


and 


-i/|  (n2-B2)  ^  ’  °>*  V-B2  (  jM|  >  uo. 


(57b) 


If  we  denote  by  I~<  the  value  of  integral  in  (54)  evaluted  for 
p  |  <  B0  and  |p|  >  BQ  respectively,  then  we  can  write 


(B2-  A1'4  ei£<d. 


-2< 


°J  [(k0o)2-b2]1/4  [b„pv-  i(b2-  m2)1/2Qv] 


(56a) 


1 


h>  “  Boe 


iit/4  f 


,  2  „2v 1/4  if>  . 

(p  -  B0)  e  d|_: 


J  [(k0p)2-M2]1/4[B0Pv+  (k2-Bo)1/2Qv] 


(58b) 


where 


13  + 


-  n2  -  n/b!” -  u2  -  p(arc  cos  -J4-  -  arc  cos  H— ) 

k0p  B0 


(59a) 


/  2  2 

f>  “  i(p*  +  'AkoP)2-  P2  *  M  arc  cos  j^)-  uln(U  ---  )  +  v^,2_  B2 


(59b) 


Setting  =  0  determines  the  stationary  points  of  the  integrals 
It  will  follow  that  the  stationary  points  will  be  given  by 


|is  =  +  B0  sin  $ 


which  are  within  the  limits  of  l2<  but  not  of  l2>  •  Hence,  using  the  method 
of  stationary  phase,  the  integral  I2  is  evaluated  to  give 


B0  cos  0ei(koP-BOcos  0  -  n/4) 


2  V k0p  Pv(4s)  -  i  cos  ^  Qv(ms) 


with 


I  Bo  sin^$  +  p4 

1  +  p 
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Comparing  (56)  and  (61)  it  will  be  noticed  that  the  contribution 
of  the  former  is  negligible  in  relation  to  the  latter,  hence  (52)  becomes 


_  -(.H-nv  pT .  cos  »Et(k°°  -  80  cos  e  '  */4> 
**30  naZD  V  itk0p  Pv(ns)  '  1  cos 


Two  partical  checks,  b  =  a  and  p  =  0,  can  be  performed  on  (63)  to 
obtain  the  well-known  result  for  a  bare  cylinder.  In  the  first  case,  b  =  a, 
one  readily  obtains 


Pv0*s)  =  0 


which  together  with  (63)  gives  the  well-known  result  [wait,  1959] 


„  y_l/Jfo  i(koP  -  AqCOs4>  -  rt/4) 
Z  If  2np  e 


obtains 


For  the  case  p  =  0,  noting  that  (BQ) (BQ)  =  i  cos  <t>,  one 

s  s 


P,,(uc)  -  i  COS  0Qv(mJ 


„U)\A 

2i 


v  s  nBo  n(i)  (Bo) 


-2  cos  » 

tta„ 


where  the  last  approximation  is  valid  if  both  A0  and  BQ  are  large  and  not 
far  apart.  Substitution  of  (66)  into  (S3)  again  yields  (65)  as  it  should, 


5.  Numerical  Results  and  Interpretation 
The  power  radiation  patterns  related  to  the  field  formula  (63)  are 
shown  in  Figs.  3  to  5 .  It  will  be  noted  that  each  pattern  exhibits  peaks  at 


certain  angles  that  appear  to  be  a  function  of  the  cylinder  radius,  Aq, 


inhoroogeneity  gradient,  p,  and  the  thickness  of  the  layer  in  wavelength  T  . 
These  patterns  resemble  in  many  respects  the  radiation  patterns  of  the  corre¬ 
sponding  planer  geometry,  homogeneous  layer  problem  considered  previously  by 
several  authors  [Tamir  and  Oliner,  1962;  Omura,  1962;  Tyras  et  al,  1964]  and 
the  curved  homogeneous  layer  considered  by  [Hasserjian,  1965]. 

The  presence  of  the  peaks  in  the  radiation  patterns  as  well  as  their 
relative  position  which  is  dependent  on  the  geometry  of  the  problem,  can  be 
explained  qualitatively  by  the  elementary  consideration  of  geometric  optics. 
Since  all  of  the  rays  radiated  by  the  source  must  satisfy  the  Snell's  law 

Je ( p)  sin  $(p)  =  Const.  (67) 

any  ray  launched  at  an  angle  o(a)  by  the  slot  will  be  deflected  while  travel¬ 
ling  through  the  layer  according  to 

D 

o(p)  arc  sin  sin  $(a)]  .  (68) 

Consequently,  the  angle  that  any  ray  will  emerge  from  the  layer,  *>(b),  will 
be  given  by 

p 

$(b)  «  arc  sin  [  sin  $(a)]  .  (69) 

It  will  follow  that  the  angle  $(b)  increases  with  increasing  radius  of  the 
cylinder,  a,  and  it  decreases  with  increasing  outer  radius  of  the  layer,  b,  and 
an  increasing  gradient  of  the  inhoroogeneity  as  measured  by  the  parameter  'p'. 
Moreover,  since  all  of  the  rays  launched  by  the  slot  are  deflected  in  the  laver 
toward  the  vertical  axis  according  to  (68),  it  will  follow  that  there  will  be 
little,  if  any,  broadside  radiation  by  this  structure.  The  pictorial  repre¬ 
sentation  of  the  foregoing,  qualitative,  analysis  is  shown  in  Fig.  6. 

The  examination  of  the  various  power  radiation  patterns  shown  in 
Figs.  3  to  5  reveals  good  agreement  with  the  qualitative  arguments  of  geometric 
optics.  In  particular,  the  presence  of  the  minor  side  lobes  in  Fig.  5  must  be 
attributed  to  the  phase  interference  effects  of  some  emerging  rays  that  occur 
under  a  certain  set  of  conditions. 


Inhomogeneous 


sheath  profile  (decreasing) 


1 
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7.  Appendix 
The  Roots  of  A,,  °  0 

From  (22)  we  have  in  this  case 


aff'ty  pv 

H^>  (B0)  ‘  Qv 

we  note 


B  -  A  *  tB 


t  a  1 


(A-l) 


(A-2) 


When  the  radii  of  the  cylinder  and  the  coating  are  large  and  not 
far  apart,  then  the  quantity  *t*  in  (A-2)  will  necessarily  be  small.  In 
this  case  we  can  use  the  approximation 

CV(A)  ~  rv(B)  -  tB  t'(B)  (A-3) 

where  £v(  )  is  any  Bessel  function.  Using  (A-3)  in  (A-l)  one  obtains  the 
following  approximate  expression 


(Bo> 


1 

(P+1)B0 


(A-4) 


Note  that  when  a  »  b  then  the  right  side  of  (A-4)  is  zero  and  the  expression 
(1)  * 

reduces  H  (B0)  »  0  which  is  an  appropriate  to  a  bare  cylinder. 

^n 

It  is  easy  to  show  using  asymptotic  form  for  the  Hankel  functions 
that  (A-4)  will  have  no  solutions  if  p  <  B0  . 

In  the  range  p  ~  BQ,  where  BQ  is  large,  one  can  use  the  transition 
region  expansion  of  the  Hankel  functions  in  terms  of  Airy  functions 
[flonl  et  al,  196 l]  to  obtain 
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<n'<V  (2  ^  A1’  t<fe)l/3«-ta/3<»»-  **»)] 

~7[s -  =  -  —  - — — “ — 7-7 - 7  •  (Ao) 

A1/  fa  \  \aol 


<V 


A1  [(|-)1/3e-i,</3(B0-  un)] 
Bo 


Substituting  (A-5)  into  (A-4)  we  write 


S,  (Bp2-  »,2)e't2,/3  _  pt!  Ai'[(^)1/3  e-U,ht0-  m)] 

(p+1)  B2/3  21/3  S  “  A1  [(|q1/3  e"I*75<V  »>] 

Bn 


(A- 6) 


The  left  side  of  (A-6)  is  small  if  ’a1  and  ’b*  are  large  and  not  too  far 
apart.  Thus  to  the  first  approximation  the  roots  will  be  determined  from 
Ai*(an)  =  0  and  consequently 


_  !/®0\l/3  in/3 

^n  ~  ‘  aA<~>  e 


(A-7) 


which  are  the  roots  appropriate  to  the  bare  cylinder. 

Since  the  first  root,  is  close  to  B0,  a  better  approximation 

to  its  value  can  be  obtained  by  expanding  the  Airy  function  and  its  derivative 
in  a  power  series.  We  note 


Ai(z) 


£l_ 

2C2 


(1  +  rr  z  -  z2  +  .  . 


) 


(A-8) 


2 

where  C2/C^  ~  0.72  [olver,  1962].  Setting  B0  -  p^  ~  2B0(B0~p^)  and  collect¬ 
ing  the  above  results,  we  obtain  a  better  approximation  for  the  first  root  as 
follows 

1/3 


,2_  £2  .  i2n/3 

“i  ~  b0  +  T7~m 

•C2l2/2  ' 2' 3 


§)y  •(P+l)et’l/3  +  2[l-(f)P+1] 


(A-9) 


i 


31 


REFERENCES 

Dettman,  J.  W.  (1962),  Mathematical  methods  in  physics  and  engineering, 
(McGraw-Hill,  Inc.,  New  York) . 

Hamm,  J.  M.,  and  G.  Tyras  (July  1965),  Further  experimental  study  of  plasma 
sheath  effects  on  antennas,  Engineering  Experiment  Station,  U.  of  Ariz., 
Rep.  No.  2,  AFCRL-65-608,  Contract  AF  19 (628) -3834 . 

Harris,  J.  H.,  Radiation  pattern  characteristics  of  line  sources  in  the 
presence  of  dielectric  layers,  Hughes  Aircraft  Co.,  Rep.  No.  P65-91, 
AFCRL-65-666  Contract  19(628) -4349 . 

Harris,  J.  H.,  and  J.  Pachares  (June  1965),  The  computation  of  cylindrical 
radiation  problems,  Hughes  Aircraft  Co.,  Rep.  No.  P65-79,  AFCRL-65-610, 
Contract  AF  19(b2B) -4349 . 

Hasserjian,  G.  (May,  1965),  Fields  of  a  curved  plasma  layer  excited  by  a 
slot,  IEEE  Trans.  Antennas  and  Propagation,  AP-13,  389-395. 

Honl,  H.,  A.  W.  Maue  and  K.  Westpfahl  (1961),  Theorie  der  Beugung,  Handbuch 
der  Physik,  vol.  XXV/1,  edited  by  S.  Flugge,  Springer -Ver lag,  p.  504. 

Liu,  Chao  Han  and  L.  Wetzel  (June,  1965),  Reflection  electromagnetic  waves 
from  inhomogeneous  plasma,  Brown  Univ.  Sc.  Rep.  No.  6,  AFCRL-65-334, 
Contr.  No.  19 (628) -2496. 

Olver,  F.  V/.  J.  (June  1962),  Bessel  functions  of  integral  order,  Handbook 
of  Mathematical  Functions,  edited  by  M.  Abramowitz  and  I,  A.  Stegun, 

NBS  App.  Math.  Series,  55. 

Omura,  M.  (December,  1962),  Radiation  pattern  of  a  slit  in  a  ground  plane 
covered  by  a  plasma  layer,  AFCRL-62-958 . 

Rusch,  W.  V.  T.  (Jan.  1964),  Radiation  from  an  axially  slotted  cylinder 

with  a  radially  inhomogeneous  plasma  coating,  Can.  J.  Phys.  42,  26-42. 

Tamir,  T.,  and  A.  A.  Oliner,  (January,  1962)  The  influence  of  complex  waves 
on  the  radiation  field  of  a  slot-excited  plaama  layer,  IRE  Trans. 
Antennas  and  Propagation  AP-10  55-65. 

Tyras,  G.,  P.  C.  Bargeliotes,  J.  M.  Hamm  and  R,  R.  Schell,  (June  1965), 

An  experimental  study  of  plasma  sheath  effects  on  antennas.  Radio 
Science  Jour.  Res.  NBS/USNC-URSI,  69D,  839-850. 

Wait,  J.  R.  (1959),  Electromagnetic  radiation  from  cylindrical  structures, 
(Pergamon  Press) . 


Unc lassifie^ 
Secunty  Classification 


DOCUMENT  CONTROL  DATA  -  R&D 

(Security  eleeei  tication  of  title  body  of  abstract  and  indexing  tnnoih  ft  on  muif  be  entered  when  the  orereJt  report  im  clee  eitiedi 
"i  ORIGINATING  AC  T I  V(  T  Y  (Corporate  tuthor,  j  2  «  R  E  PO»  T  SE  C  uB,  T  y  C  l*55Tic*tion 

Engineering  Experiment  Station  j  Unclassified 

The  University  of  Arizona  f  2b  GROUP 

Tucson,  Arizona  j 

3  REPORT  TITLE 

Field  of  an  Axially  Slotted  Circular  Cylinder  Clad 
with  an  Inhomogeneous  Dielectric 


4  DESCRIPTIVE  NOTES  ( fype  oi  report  and  tnctueive  date*) 

Scientific  Report  (Interim)  August,  1965  -  January,  1966 

5  AU  THORfSj  (Lett  name,  ft  ret  nemo,  initial) 

Tyras,  George 


6  REPO  RT  OATE 

February  1966 

7a  TOTAL  no  of  pacej 

38 

7b  NO  OF  REFS 

13 

8a  CONTRACT  OR  GRANT  NO 

AF  19(628) -3834 

b  project  no  ana  Task  No.  4642,02 

Sa  ORIGINATOR'S  REPORT  NUMBERS; 

No.  3 

c  DOD  element  62405  304 

9b  other  report  no(S)  (A  ny  other  numbere  that  mey  be  seel gned 

Chi*  report) 

a  DOD  subelement  674642 

. 

AFCRL-66-203 

10  availability/limitation  notices 

Distribution  of  this  document  is  unlimited. 

1<  SUPPL EMENTARY  NOTES 

... 

12  SPONSORING  MILITARY  ACTIVITY 

Hq.  AFCRL,  OAR  (CRD) 

United  States  Air  Force 

L.  G.  Hanscom  Field,  Bedford,  Mass. 

13  ABSTRACT 

Field  expressions  appropriate  to  small  and  large  cylinders  are  found 
when  the  cylinder  is  coated  with  a  layer  of  inhomogeneous  dielectric  of  the  form 
€ (p)  «»  (p/b)2p,  a  <  p  <  b,  where  ’a'  is  the  radius  of  the  conducting  cylinder, 

* b *  is  the  outer  radius  of  the  dielectric  coating  and  * p *  is  an  arbitrary 
parameter.  (U) 

In  the  case  of  large  cylinders  coated  with  plasma,  the  radiation 
patterns  are  plotted  for  various  combinations  of  the  cylinder  radius,  thickness 
of  the  coating  and  the  inhomogeneity  gradient.  The  radiation  patterns  are  in 
good  agreement  with  qualitative  arguments  of  geometric  optics. (V) 


form 

1  JAN  C4 


1473 


0101 -807-6800 


Unclassified 

Security  Classification 


Y 

i 


Unclassified _ 

Security  Classification 


14 


KEY  WORDS 


Slot  Antenna 
Plasma  Sheath 
Radiation  Pattern 
Inhomogeneous  Dielectric 
Inhomogeneous  Plasma 


LINK  A 


ROL  E 


LINK  0 


LINK  C 


ROLE 


T 


INSTRUCTIONS 

t.  ORIGINATING  ACTIVITY  Enter  the  name  and  address 
of  the  contractor,  subcontractor,  grantee.  Department  of  De¬ 
fense  activity  or  other  organization  < co-porate  author)  issuing 
the  report, 

2a.  REPORT  SECURITY  CLASSIFICATION:  Enter  the  over¬ 
all  security  classification  of  the  report.  Indicate  whether 
"Restricted  Data"  is  included.  Marking  is  to  be  in  accord¬ 
ance  with  appropriate  security  regulations. 

2b.  GROUP:  Automatic  downgrading  is  specified  in  DoD  Di¬ 
rective  5200.10  and  Armed  Forces  Industrial  Manual.  Enter 
the  group  number.  Also,  when  applicable,  show  that  optional 
markings  have  been  used  for  Group  3  and  Group  4  as  author¬ 
ized. 


imposed  by  security  classification,  _s.ng  standard  statements 
such  as: 


3.  REPORT  TITLE:  Enter  the  complete  report  title  tn  all 
capital  letters.  Titles  in  all  cases  should  be  unclassified. 

If  a  meaningful  title  cannot  be  selected  without  classifica¬ 
tion,  show  title  classification  m  all  capitals  in  parenthesis 
immediately  following  the  title. 

4.  DESCRIPTIVE  NOTES  If  appropriate,  enter  the  type  of 
report,  e.g. ,  interim,  progress,  summary,  annual,  or  final. 

Give  the  inclusive  dates  when  a  specific  reporting  period  is 
covered. 

5.  AUTHOR(S).  Enter  the  name(s)  of  authorfs)  as  shown  :n 
or  in  the  report.  Entei  last  name,  first  name,  middle  initial. 

If  military,  show  rank  and  branch  of  service.  The  name  of 
the  principal  author  is  an  absolute  minimum  requirement. 

6.  REPORT  DATE.  Enter  the  date  of  the  report  as  day, 
month,  year,  or  month,  year.  If  more  than  one  date  appears 
on  the  report,  use  date  of  publication. 

7 a  TOTAL  NUMBER  OF  PAGES  The  total  page  count 
should  follow  normal  pagination  proi  edures.  1.  e. ,  enter  the 
number  of  pages  containing  information. 

76  NUMBER  OF  REFERENCES  Enter  the  total  number  of 
references  cited  in  the  report. 

So  CONTRACT  OR  GRANT  NUMBER  If  appropriate,  enter 
the  applicable  number  of  the  contract  or  grant  under  which 
the  report  was  written. 

86.  8c.  h  8 d  PROJECT  NUMBER  Enter  'ne  appropriate 
military  department  identification,  such  as  project  number, 
subproject  number,  Bystem  numbers,  task  number,  etc. 

9 a  ORIGINATOR’S  REPORT  NUMBER(S)  Enter  the  offi¬ 
cial  report  number  by  which  the  document  will  be  identified 
and  controlled  by  the  originating  activity.  This  number  must 
be  unique  to  this  report. 

96  OTHER  REPORT  NUMBER(S)'  If  the  report  has  been 
assigned  any  other  report  numbers  f either  by  the  origin otor 
or  by  the  sponsor),  also  enter  this  number(s). 

10.  AVAILABILITY  ^LIMITATION  NOTICES'  Enter  any  lim¬ 
itations  on  further  dissemination  of  the  report,  other  than  those 


(1) 


(2) 


(3) 


(4) 


(5) 


"Qualified  requesters  may  obtain  copies  of  this 
report  from  DDC.  ” 

"Foreign  announcement  and  dissemination  of  this 
report  by  DDC  is  not  authorized.” 

"U.  S.  Government  agencies  may  obtain  copies  of 
this  report  directly  from  DDC.  Other  qualified  DDC 
users  sh.il  1  request  through 


"U.  S.  military  agencies  may  obtain  copies  of  this 
report  directly  from  DDC  Other  qualified  users 
shall  request  through 


"All  distribution  of  this  report  is  controlled.  Qual¬ 
ified  DDC  users  shall  request  through 


If  the  report  has  been  furnished  to  the  Office  ol  Technical 
Services,  Department  of  Commerce,  for  sale  to  the  public,  indi¬ 
cate  this  fact  and  enter  the  price,  if  known. 

11.  SUPPLEMENTARY  NOTES  Use  for  additional  explana¬ 
tory  notes- 

12.  SPONSORING  MILITARY  ACTIVITY.  Enter  the  name  of 
the  departmental  project  office  or  laboratory  sponsoring  Spay¬ 
ing  ler)  the  research  and  development.  Include  address. 

13  ABSTRACT.  Enter  an  abstract  giving  a  brief  and  factual 
summary  of  the  document  indicative  of  the  report,  even  though 
it  may  also  appear  elsewhere  in  the  body  cf  the  technical  re¬ 
port  If  additional  space  is  required,  a  continuation  sheet  shall 
be  attached 

It  is  highly  desirable  that  the  abstract  of  c  lassifiea  reports 
be  unclassified  Each  paragraph  of  the  abstract  shall  end  with 
an  indication  of  the  mihtarv  security  classification  uf  the  in 
formation  in  the  paragraph  represented  as  (TSj  -'Si  'C ,  nr  (U) 

There  is  no  limitation  on  the  length  of  the  abstract  How¬ 
ever,  the  suggested  length  is  from  1  50  !  >  226  words 

14  KEY  WORDS  Key  words  are  technically  meaningful  terms 
or  short  phrases  that  characterize  a  report  and  may  be  used  as 
index  entTies  for  cataloging  the  report  Key  words  mu?>!  be 
selected  so  that  no  security  classification  is  required  Identi 
fiers,  su<  h  as  equipment  model  designation,  trade  name,  military 
project  code  name,  geographic  location,  may  be  used  as  key 
words  but  will  be  followed  by  an  indu  anon  of  technical  con¬ 
text  The  assignment  of  links,  roles,  and  weights  is  optional 


Unclassified 


Security  Classification 


re-  "T" 


“f 


